IMM-NYU  324 
APRIL  1964 


NEW  YORK  UNIVERSITY 
COURANT  INSTITUTE  OF 
MATHEMATICAL  SCIENCES 


Buckling 


EDWARD  L.  REISS 


raEPARED  UNDER 

<^ANT  NO.  DA*AR0-(D)-5M240344 
WITH  THE 

U.S.  ARMY  RESEARCH  CH^FICE 
AND 

CXJNIRACr  2«>.  NONR-285(42) 
WriHtHfi 


IMM-NYU  32k 
April,  1964 


New  York  University 

Courent  Institute  of  Matheroatical  Sciences 


BIFURCATION  BUCKLING  OF  SPHLRICAL  CAPS 

Edward  L.  Re las 


This  report  represents  results  obtained  at  tbe  Courant 
Institute  of  Hathenatlcal  Sciences,  New  York  University, 
under  the  sponsorship  of  the  U»S*  Army  Research  Office, 
Grant  No.  DA-ABO<-(D) -31  **124-0 344  ®nd  the  Office  of  Naval 
Reaearch,  Contract  No.  Nonr-265(42} •  Reproduction  in 
whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  GoTemwent. 


1 


Abstract 

A  nonlinear  boundary  value  problem  is  considered  for 
the  axisyrametric  buckling  of  thin  spherical  shells  subjected 
to  uniform  external  pressure.  The  unifoxnaly  compressed 
spherical  state  is  a  solution  of  this  problem  for  all  values 
cf  the  pressure.  We  prove,  using  Poincar6»s  method,  that 
for  pressures  sufficiently  near  each  simple  eigenvalue  of 
the  linearized  shell  buckling  theory,  there  is  another 
(buckled)  solution  of  the  nonlinear  problem.  A  convergent 
perturbation  expansion  is  used  to  analyze  the  buckled  solutions 
near  the  eigenvalues.  For  a  limited  range  of  caps,  we  also 
prove  that  one  or  three  buckled  solutions  bifurcate  from  the 
multiple  (double)  eigenvalues  depending  on  their  order.  The 
existence  of  a  ’’lowest”  intemedlate  buckling  is  established 
and  precise  upper  and  lower  bounds  are  given  on  its  wignitude. 

Introduction 

The  surface  of  a  thin  elastic  spherical  cap  is  subjected 
to  a  uniform  pressure,  p,  which  is  directed  towards  the  cap ^3 
center  of  curvature*  It  is  a  well  known  experimental  result, 
see  e*g«  [l] ,  that  as  the  pressure  increases  from  zero  the 
cap  defoms  only  slightly  from  the  spherical  shape  until  a 
critical  pressure  p  «  is  reached.  Then  the  cap  suddenly 
jumps,  with  relatively  large  deflections,  into  a  non-spherlcal 
shape  which  we  call  the  buckled  state*  The  fundamental  problem 
of  shell  buckling  is  to  determine  the  mechanism  which  initiates 
the  jumping  and  to  obtain  estimates  of  p^. 
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Previous  investigators  «  have  assumed,  as  on  approxima¬ 
tion  to  experimental  conditions,  that  the  cap's  edge  is  rigidly 
clamped,  i.e»  the  displacement  and  chrnge  of  slope  are  zero# 

In  this  paper  edge  conditions  are  considered  for  which  the 
spherical  shape  (the  unbuckled  solution)  is  a  pos'!ible  solution 
of  the  nonlinear  problem  for  all  pressures*  Vis  refer  to  these 
as  oifurcation  buckling  problems  since  other  (bu-:kled)  solu¬ 
tions  of  the  nonl'. near  problem  nay  branch  from  the  unbuckled 

A 

solution*  For  example,  a  bifurcation  buckling  problem,  which 
we  call  ^roblen  B,  is  obtained  if  the  follov.’ing  conditions  are 
specified  on  tho  edge  of  the  cap:  no  rotation  (clamped); 
zero  transverse  shear  force,  i*e«,  the  edge  is  free  to  move 
normal  to  the  spherical  surface;  and  the  meridional  membrane 
stress  is  prescribed  so  that  it  is  in  equilibrium  with  the 
applied  surface  pressure.  Other  bifurcation  problems  are 
obtained,  for  exaripie,  by  replacing  the  condition  on  the 
meridional  stress  in  Problem  B  with  a  corresponding  one  on  t^e 
meridional  membrane  displacement,  or  by  permitting  the  cap  to 
frecAj  rotate  instead  of  clamping  it.  The  bifurcation  problems 
are  precisely  formulated  in  Section  2#  In  Section?  3,  4  sjid  5 
only  Problem  3  is  considered.  However  analogous  results  can 
be  established  for  other  bifurcation  problems.  Some  of  these 
are  contained  in  the  fins!  section  of  the  paper.  We  consider 


^  See  e.g,  [l»7]  and  references  contained  therein. 

^  For  the  rigidly  clamped  cap,  which  wo  refer  to  as  relaxation 
buckling,  the  spherical  shape  is  a  solution  if  and  only  if 

p  *  0. 
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only  axiaymmetric  deformations  of  the  cap. 

Bifurcation  buckling  problems  for  spherical  cap^  were 
first  considered  in  an  approximate  form  in  [8] ,  An  equivalent 
of  Problem  B  was  treated  in  [5]  where  the  lineari-^^d  buckling 
theory  was  partially  analyzed  and  approximate  solutions  of 
the  nonlinear  problem  were  obtained#  The  linear  buckling 
theory  had  been  previously  discussed.* 

The  precise  knowledge  of  the  unbuckled  state  for  the 
bifurcation  problems,  permits  us  to  rigorously  establish 
certain  properties  of  the  solution.  In  Section  3  we  prove 
that  for  all  P  sufficiently  near  each  simple  eigenvalue  of 
the  linearized  shell  buckling  theoryt  a  solution  of  the  non¬ 
linear  problem  exists#  Only  one  solution  of  ttu  nonlinear 
problem  branches  from  each  simple  eigenvalue#  A  related  result 
is  established  in  Section  U.  for  t*-s  other  eigenvalue^,  which 
are  all  double,  and  for  a  limited  range  of  caps#  However,  we 
find  the  surprising  result  that  one  or  three  solutions  of  the 
nonlinear  problem  bifurcate  from  the  double  eigenvalues 
depending  on  their  order#  The  conjectured  load-de format ion 
curves  for  simple  eigenvalues  ere  sketched  in  Pig.  1.  A 
perturbation  expansion,  which  is  valid  near  each  eigenvalue, 
is  used  to  prove  that  the  curves  in  Pig.  1  have  the  form  of 
the  solid  portion# 

Our  method  of  analysis  is  the  Poincar^  bifurcation 
theory  u^ed  to  prove  the  existence  of  periodic  solutions  of 

*  '  ' .  '  '  .  II" 

Privately  coxanunicated  to  the  author  by  W.  Squires,  1958. 

f 

'  Here  P  is  a  diaenaionleaa  parasaeter  proportional  to  p,  see 
Eq#  (2»2a}  below# 
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initial  value  problems.  Previous  applications  of  this  tech¬ 
nique  to  boundary  value  problems  are  given  in  [9  f  lO]  where 
related  buckling  problems  for  columns  and  circular  plates  a’-e 
studied.  In  Section  3»  where  t  e  simple  eigenvalues  are 
studied,  the  procedures  employed  are  closely  related  to  those 
given  in  [lO] •  Modifications  are  made  in  Section  U  to  inves¬ 
tigate  the  double  eigenvalue  case# 

Friedrichs  [ll]  proposed  an  energy  mechanism  and  intro¬ 
duced  the  concept  of  an  intermediate  buckling  load  to  explain 
the  experimentally  obse^^ved  juraping  of  con!)lete  spheres  from 
an  unbuckled  to  a  buckled  state.  Modifications  and  extensions 
of  these  ideas  were  subsequently  proposed  in  [l2].  As  applied 
to  the  bifurcation  buckling  of  caps,  see  e.g.  Pig.  1,  sev- 
eral  intermediate  buckling  loads  '  corresponding  to  differ¬ 
ent  branches  of  the  solution  Diay  exist.  For  a  given  branch, 
say  bifujrcating  from  the  eigenvalue  P  ,  pA*^^  is  defined  as 
a  load  in  the  interval  <P  such  that  for  all  P  In 

p^n)  <  p  <  the  unbuckled  state  has  less  potential  energy 

than  the  corresponding  buckled  state  and  conversely  for  P  in 

<  P  <  P  •  In  Section  5  we  establish,  assuring  that  the 
n  n 

potential  energy  has  a  riinimun  for  every  finite  P  and  p  >  0, 


the  existence  of  a  lowest  intermediate  buckling  load  see 
Pig.  1.  It  is  shown  to  be  bounded  from  belw  by  the  lowest 
buckling  load  of  an  "equi  lent"  flat  circular  plate  buckling 
problem  and  bounded  from  above  by  the  lowest  eif-*nvalue,  of 
the  llneariaed  shell  buckling  theory.  More  accurate  upper 
bounds f  which  are  considerably  less  than  P,  are  rigorously 


5 


obtained  in  Section  5  by  a  minimization  procedure.  These 
bounds  bracket  for  a  limited  range  of  caps.  It  seems 
likely  that  with  suitable  modifications  corresponding  esti¬ 
mates  of  intermediate  buckling  loads  can  be  obtained  for  the 
unsyrametric  buckling  of  spherical  caps  and  other  shell  bifur¬ 
cation  buckling  problems.  Some  of  the  results  of  Section  5 
are  closely  related  to  ones  previously  announced  by 
Vorovich  []l3]  • 

5,  Formulation  of  the  Boundary  Value  Problem 

The  elastic  spherical  cap  is  of  thickness  2h  and  radius 
R  and  has  a  small  angle  of  opening  Z/\  ;  see  Pig,  2  for  the 
shell  geometry.  We  consider  the  axisymmetric  deformations  of 
the  cap  that  result  frcan  a  uniform  and  inwardly  directed  pres¬ 
sure  p.  The  non-vanishing  mid-surface  displacements,  u  and  w 
which  are  in  the  meridional  and  normal  directions  to  the  sheila 
mid-surface  are  therefore  functions  only  of  the  polar  angle  0. 
Both  w  and  p  are  counted  positive  when  directed  towards  the 
center  of  curvature.  Nonlinear  differential  equations  which 
describe  the  small  finite  axisymmetric  deformations  of  thin 
spherical  caps  have  been  derived  by  several  authors,  e.g, 
[1,2,14«*16],  These  equations  may  be  written  as, 

(2.U)  at'.x)  ♦xf{»)  «  f«(x)[rtx)  +1]  , 

(2.1b)  Og(x)  •  -f»[f^(x)  ♦  2f (x)J  . 
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Here  we  have  used  the  dimensionless  variables: 

X  =  9//,  ,  f(x)  =  .  p  =  (A^/C)(R/h)  . 

(2.2a) 

F  =  (R/h)^(p/22;C)  ,  -K  ^  ?P  ,  ~  2/3(1  -v^)  , 


and  the  linear  differential  operator  G  which  is  defined  by 

G^(x)  =  x"^[x^^ ’ (x)] *  , 


where  a  prime  indicates  differentiation  with  respect  to  x. 

Here  E  is  Young’s  modulus  and  v  is  Poisson’s  ratio, 

The  "excess"  stress  function  g(x)  in  (2.1)  is  defined 
such  that  the  meridional  and  c ire ujtife rent ial  membrane  stresses 
Oq  and  0 

are  given  by 


^  and  the  corresponding  dimensionless  stresses  ^  and 


[[^(x)  =  (R/h)  (2/EC)Oq(9)  =  g(x)  -P  , 

(2.2b) 

J^(x)  5  (h/h) (2/EC)o^(©)  =  [xg(x)  -xP]  ’  . 

The  outer  surface  bending  stresses  and  c|,  and  the  corre¬ 
sponding  dimensionless  stresses  and  are  given  in  terms 
of  f(x)  by 


=  (R/h)(2/3EC^)o°(9)  =  xf’(x)  Ml*  v)f(x)  , 

(2.2c) 

^(x)  5  (R/h)(2/32c2)o|(9)  =  vxf'(x)  ^  (!<►  v)f(x)  . 


The  independent  variable  f(x)  is  related  to  the  slope  of 
the  deformed  middle  surface  of  the  cap  with  respect  to  the 
initial  spherical  shape.  Thus  if  f(x)  =  0  for  a  given  defor¬ 
mation  the  deformed  middle  surface  is  also  spherical.  We 

^  *  -  -  -  -  — 

The  independent  variables  a{x)  and  y(x)  employed  in  previous 

papers  [2,3]  are  related  to  the  present  variables  by. 
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refer  to  P,  which  is  defined  in  (2.2a),  as  the  f^eonietr ical 

I 

parameter  and  P  and  X  as  either  pressure  or  loading  parameters. 

To  complete  the  fo.’nnulation  conditions  at  the  center 
X  =  0  and  the  edge  x  =  1  are  required.  From  the  aymnetry  of 
the  deformation  and  the  regularity  of  the  membrane  and  oend- 
ing  stresses  a^-  the  origin  we  obtain  with  aid  of  (2.1)  that, 

(2«3)  f'(C)  =  g’(0)  =  C  . 

The  edge  of  the  cap  is  restrained  from  rotating  so  that, 

(2.4a)  1*(1)  =  0  . 

In  addition,  we  assume  that  at  x  =  1  the  transverse  shear 
force  vanishes  and  5^(1)  is  specified  so  t  .at  the  cap  is  in 
equilibrium  with  the  applied  pressure.  This  yields  the  bound¬ 
ary  condition 

(2.4b)  g(l)  =  C  . 

The  bifurcation  buckling  problem.  Problem  B,  is  defined 
as  the  boundary  value  problem  consisting  of  the  differential 
equations  (2.1)  and  the  boundary  conditions  (2.3)  and  (2,4) • 
Other  boundary  conditions  can  be  specified  at  x  =  1 
to  yielc*  bifurcation  buckling  problems  for  clamped  caps,  e.g. 


(2.5)  f(l)  «  0 


* 


g'  (1)  ♦  (1  -  v;g(i)  =  c  . 
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These  conditions  imply  that  on  the  edge  the  transverse  shear 
force  vanishes  and  the  meridional  displacement  u  (or  equiva¬ 
lently  the  horizontal  displacement)  is  proportional  to  P  The 
proportionality  constant  is  determined  so  that  the  spherical 
form  is  a  possible  solution  for  all  P, 

3.  The  Existence  of  Buckled  Solutions. 

A  solution  of  Problem  B  that  is  valid  for  all  finite  P 
and  p  is  the  unbuckled  solution,  f(x)  =  g(x)  =  0.  This  corre¬ 
sponds  to  a  state  of  uniform  compression  in  which 

Z©  (==)  =  -p . 

and  the  deformed  middle  surface  remains  spherical. 

The  existence  of  buckled  states  will  now  bo  established 
using  Poincare's  method.  Specifically,  we  prove  in  this  sec¬ 
tion  that  for  every  0  >  Q  and  for  every  positive  integer  n 
there  exists  a  buckled  solution  when  P  is  in  a  sufficiently 
small  and  full  interval  about  the  n-th  simple  eigenvalue, 
p^,  of  the  linearized  shell  buckling  theory.  This  is  in  con¬ 
trast  to  the  buckled  circular  plate  [lOj  where  the  load  must 
be  slightly  greater  than  each  eigenvalue  of  the  corresponding 
linear  buckling  theory.  V/e  also  show  that  for  P  sufficiently 
near  P^^  the  solution  has  n-1  internal  nodes.  In  the  follow¬ 
ing  section  the  solutions  of  Problem  B  near  the  other  eigen¬ 
values  of  the  linearized  theory,  which  are  all  double,  are 
studied. 
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To  employ  Poincare’s  method  parameters  c  and  ‘and  new 
independent  variables  y(x}  and  z(x)  are  defined  by, 

(3,1a)  e  =  lim  [xf  (x)]  *  ,  6=  lim  [xg(x)J  ‘  , 

X  — »  0  X  0 

(3.1b)  y(x)  =  e“^xf(x)  ,  z(x)  5  e“^xg(x)  , 

The  differential  equations  (2.1)  of  Problem  B  are  then  given 
by, 

(3,2a)  hy(x)  +Ppy(x)  =  rz (x)  [l  +  £5  ;x)/x]  , 

(3.2b)  Hz(x)  =  -u[ey^'x)/x  +  2y(x)]  , 

where  h  is  the  linear  differential  operator, 

(3.2c)  =  ,x’’^[x^(x)]  ’  >  ^  . 

The  initial  value  problem,  ,  is  defined  by  the 
differential  equations  (3.2)  and  the  initial  conditions: 

(3.3a)  y(0)  =  0  ,  y’(0)  =  1  , 

(3.3b)  ziO)  =  0  ,  2'(0)  = 

Per  all  finite  and  r  >  0,  Problem  *.  has  a  unique 

solution  in  the  interval  0  ^  1  which  is  analytic  in 
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P,  t  and  S*  This  result  can  be  proved  by  toe  same  methods 

used  in  [lo]  to  establish  an  analogous  result  for  the  buckling 

of  circular  plates  and  hence  the  proof  is  not  given  here*  For 

fixed  p  >  0  values  of  the  pa^-ameters  P,  e  and  are  sought 

such  that  the  solution,  fy(x;P, c, 6)  #  2(3t;P,  e,£5)j  of  satisfies 
the  boundary  conditions 

O.li)  y(l;P»s,j)  =  z(ljP,e,>)  =  o  . 

If  such  a  choice  of  parameters  is  possible  it  then  follows 
from  (3.1b)  that  f(x)  =  and  g(x)  = 

Is  a  solution  of  5« 

A  solution  of  which  satisfies  the  boundary  conditions 
(3,14.)  la  obtained  by  choosixig  the  special  parameter  value, 
e  =  0  in  (3»2)«  The  resulting  initial  value  problem  for  the 
functions  y  =  y(x;?,0,-^),  a  =  z(xjp,0,  i)  has  solutions  which 
satisfy  (3»l4-)  ii*  and  only  if 

(3.S«)  P(f)  =  Fn(r)  5  2-/  I*  n^e  2.  2/2 
(3.5b)  ^(r)  =*„(.’)  5  2  '/;  2 

where  is  the  n-th  zero  of  the  Bessel  function  J^(x).  The 
quantities  ^he  eigenvalues  (or  buckling  loads)  of 

the  linearized  buckling  theory  and  they  are  simple  eigenvalues 

Fnm  ^  ^  » 

'  Din  n  mn  m  mn 


n  ”  1  f  2  j 
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(3.6) 


^  n-1/2 


mn 


£J 

m  n 


m 


►  A  •  *3  ^  n 


The  solutions  are  the  eigenfimctions  (or  buckling  inodes)  of 

it 

the  linearized  theory  and  for  simple  eigenvalues  are  given  by 


(3.7) 


I 

I 


f  n  a  1,2, 


•  • » 


5  2  /..h[(0) 


Tb«  corresponding  solutions  of  Problem  B  reduce  to  the  unbuck¬ 
led  solution* 

Thus  for  each  finite  p  >  0  there  are  a  denuraerably 
infinite  number  of  ro'-ts  “  “  1,2,...  of  (3*i4.). 

According  to  the  implicit  function  there  art?  other  roots  of 
(3.U)  near  appropriate  Jacobian 


(3.8)  J  s  *1"^  (Dyi”^  (1)  -  (Dy!*^^  (l)  5^  o 


n  "  1,2, #•< 


Here  we  have  used  the  notation 


For  simplicity,  the  explicit  dependence  of  the  solutions  on 
p  is  suppressed  here  and  in  the  remainder  of  the  section. 


! 
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(3.9) 


y(x;P,e,^) 


etc* 


P=P 


n 


es=0 

6=^ 


n 


1,2, 


•  •  »  • 


Since  y  and  z  are  analytic  in  P  and  <3,  the  initial  value 

problems  (variational  problems)  satisfied  by  jp,  Zp 

and  2<-  can  be  obtained  by  formal  differentiation  of  with 
b 

respect  to  F  and  o*  Thus  we  obtain  from  (3*2)  and  {3*3) 
using  (3*7)  and  (3*9), 


(3.10a)  Hy<"*+P^py<"' , 
(3.10b)  Hz*"'  =  . 

(3.11a)  Hy^"’ +  , 

(3.11b)  =  -2i'y^”*  , 

b  ' 

The  explicit  solution  of  (3*11) 
at  X  =  1,  the  relations 


y^"’(0)  =  y'"''(o)=o,’ 
z^"'(0) 

.rFa,2,*« 

y^^^O)  =y^J^'  (C)=0,| 

i 

z^.^^(O)  s0,z^’^(0)»li 
6  / 

(not  presented)  yields. 


'  ^  0  it  P  ^  P 
^  mn 

(3-12)  *2^’^y!^^  (1)  •>  ,  Bi,n=l,2,...,  a^. 

O 

*'inn 


=  0  if  P  = 


Therefor©  it  folloKs  from  {3*8)  and  (3.12)  that 
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(3.13)  J  =  y'"^l)|^^*^y^"^(l)  -2|"’(1)  . 

HeRC©f  J  «  0  if  P  =  and  P  -  P  and  nothing  can  be  concluded 

I  '  mn  n 

concerning  the  existence  of  neighboring  roots  to 

*  This  case  is  considered  in  the  following 

section*  To  prove  that  J  0  if  P  -  ,  P  P_  we  assume 

ton  n 

the  contrary,  i«e.  let, 

(3.14)  yp"*(l!  =  =  a 

vbere  a  Is  an  arbitrary  real  constant.  Then.lyl'’^  (i)  .zl"' (x)] 
is  a  solution  of  (3*10)  that  satisfies  the  boundary  conditions 
{3*lii)«  A  .solution  of  this  boundary  value  problem  exists 
if  and  only  if  the  inhor.ogsneous  terms  are  orthogonal  to  every 
solution  of  the  hoMogsneous  adjoint  problem.  This  leads  to 
th©  condition  that  for  all  a 

1 

'0 

which  Is  impossible  and  hence  J  ^  0. 

Hence,  by  the  implicit  function  theorem  and  the  analytic 
city  of  y  and  s  in  the  parameters,  (3.4)  Pc  uniquely 
solved  for  P  and  6  as  analytic  functions  of  e  In  some  suf¬ 
ficiently  ssiall  neighborhood  of  each  root  [P^»0,P'  ]  ,  c/ 

The  solutions  are  the  analytic  functions 


Ik 


(3»l5)  P  *  ^  *^n*  ^  ~  1»2,«»»  f 

which  satisfy  the  conditions 


'’(n)<°)  =  Pn 


For  sufficiently  small  e 


(3.16) 


g(x;c)  =  ex“^2^(x;6)  5  ex*'^z(x;P^^j  (e)  ,e,  (e)) 


-1. 


are  solutions  of  B. 

The  solutions  of  B  near  c  =  C  are  now  considered.  Since 
p^^j(c),  y^(xje)  and  2^(x;p'»  are  analytic  functions 

there,  they  have  convergent  expansions  in  sorae  interval 
hi  <  e®  given  by: 


„)U)  - 


^(n) 


(3.17) 


(x)^'^yj'^^  (x)  e^,  Zj^(x;e)  «  (x)^^  ^2^^Cx)c^ 


The  expansion  coefficients  are  detemined  in  the  usual  w^y  by 
substituting  (3.17)  into  (3.2)  and  (3.U)  and  equating  the 
coefficients  of  each  power  of  e.  Bach  of  the  resulting 
system  of  linear  boundary  value  problens  for  the  coefficients 

[y^*^,  has  a  solution  if  and  only  if  the  appropriate 


orttiogonality  condition  is  satisfied*  For  the  ^irst  two  of 
these  problems  with  i  =s  1,2  the  orthogonality  conditions 
yield 


Thus  it  follod’s  from  the  first  of  (3»17)  and  (3*l8a}  that  for 
each  n  solutions  of  B  must  exist  for  all  P  in  some  small 
intex^l  about  P^. 

For  sufficiently  small  c  we  have  from  (3*17) 


Benee  in  the  interval  0  <  x  <  1,  y^  has  the  same  number  of 

seros  as  or  y  has  n-1  simple  teros* 

o 

Approximations  of  the  lower  buckling  load, 
for  each  branch  may  be  obtained  by  truncating  the  series  in 
the  first  of  (3«17}  and  from  the  condition  that 

— — -0  •»  the  lower  buckling  load.  These  approximations 
may  be  valid  only  if  the  expansions  (3*17)  converge  for  suf¬ 
ficiently  large  c*  For  example  if  2  terras  are  retained  in  the 
series  so  that 
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then 


Pj^(  e)  *►  P^’’^  e^ 


p|n)2 

•’(n)L“*’n'y,W  • 


4*  The  Existence  of  Buckled  Solutions  Near  Double  Elgenvaluea 

We  now  consider  P=  and  the  ^oubie  eigenvalues, 

ton 


(U.l)  P  = 


P  (  '  )  = 

m  :nn 


n 


'^mn> 


m,n  =  1,2,...,  XL  /  n  • 


The  solutions  of  the  initial  value  problem  (3*2),  (3*3)  with 

e  ^  0,  r  =  "  P  -  P  ('-  ),  that  satisfy  the  boundar’r  con- 

'  mn  m  mn 

ditions  (3*i;)  are  the  eigenfunctions'^ 


{U*2a) 


where  ^  is  an  arbitrary  real  nunber  and 


(U.2b) 


mn 
ol/2  » 


a  m  - - 3 - 


-1] 


Kere,  end  in  all  subsequent  eroationa  in  this  section, 
t  e  Indicias  m,n  ”  1,2,...  ,  n  n,  and  we  shall  hereafter 
omit  explicit  reference  to  this. 
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The  double  eigenvalue  case  may  be  cf  eeoeeial  physical 
significance.  Consider  the  variation  with  p  of  ?("')  which  is 
defined  by 

(i;»3)  P(c)  5  min  Pjj{p)  * 

n 

TkMjre  are  local  ninirauias  of  P(P)  at  P  =  C  =  2^^^^’^, 

“•  n 

n  e  1,2,...  .  At  >'' «  is  discontinuous  or  there 

are  "peaks^  in  the  curve  of  P(p)  at  the  double  eigenvalues. 

A  corresponding  peaking  behavior  was  previously  noted  [2,17j[ 

in  the  experinental  buckling  loads  obtained  by  Kaplan  and 

Pung*[l].  Although  the  experimental  boundary  conditions 

probably  differ  from  those  of  Problem  3,  this  phencxaenon  may 

be  related  to  the  occurrence  of  a  double  eigenvalue. 

Since  the  Jacobian  in  (3*6)  vanishes  it  P  ^  o  and 

mn 

pap,  the  procedure  given  in  Section  3  must  be  modified  to 
investigate  this  case.  The  modifications  are  essentially 
extensions  of  the  methods  used  in  the  Poincar^  theory  when 
the  variational  problem  possesses  periodic  solutions,  see 
[l8] •  A  new  parenster  ^  is  introduced  and  modi/ led  initial 

value  probleau,  for  the  functions  y _ (x;P,e,5,k)  and 

t^(x;P,c,  ,<)  are  defined  by 


*  For  example,  the  experlsiental  peaks  occur  at 


20,  $6  and 
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Hy  ♦P.'’  y  z  s  to  x'^y  z  , 

^mn  irai‘'mn  ran  mn  inn  •^mn  ran  ■'  ’ 


(U.4) 


hZ_  *  2C>  y„  =  -^Knn^’^yfn  ' 

ran  ran  ran  ran  •'ran 


y««  =  0,  y„„  =  1,  z^^  =  0,  z^^  =  6,  for  x  =  0  . 
•'ran  •'ran  ran  ran 


We  can  show,  as  in  the  preceding  section  that  the  unique 

solution  of  ^  Is  an  analytic  function  of  the  pararaeters 
mn 

P,  e,  S  and  K.  The  parameter  values  are  to  be  determined 

such  that  the  solutions  of  satisfy  the  boundary  conditions, 

ran 

=  2^(ljF,e,5,x)  =  0  . 

For  the  special  parameter  values  e  =  X  =  0  and  P  =  P  =  P 

ran 

the  solutions  of are 

satisfy  (4.5) 

for  all  bounded  and  real  •  The  implicit  function  theorOM 

ran 

is  used  to  determine  other  roots  of  (4»5)  near  each  root 

*  Proceeding  as  in  Section  3  and  using 

the  appropriate  variational  problems  we  can  show  that  the 

Jacobian  of  (4,5}  with  respect  to  P  and  ^  and  evaluated  at 

each  root  does  not  vanish  for  all  Thua  there 

ran  inn  tbo 

are  solutions  of  (4*5}, 


P 


ran 


(4.6) 


f 
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which  «re  analytic  functions  of  c  and  6  in  some  sufficiently 

small  neighborhood  of  each  root  [?  (r’  ) ,  for  which 

ic  mn  Pin 

^mn  ^  ^^mn'  ^rm  aatlsfy  the  conditions 

Hence  the  modified  problems  (l4-«U)»  (U*5)  have  solutions 

*(iBn)  which  are  analytic 

in  e  and  6  for  sufficiently  sisall  |ej  and 

Solutions  of  the  modified  problems  are  solutions  of 
the  original  problem  (3*2),  (3.3)  and  (3*U)  if  6  and  e  are 
chosen  so  that  the  bifurcation  eouations  [l6j 

(U.«  8 )  t  ( t ,  ^)  ~  0  , 

WPt  satisfied*  Thus  if  (U*8)  have  solutions  6  ^  ^  (c)  v#ilch 

mn 

satisfy  the  ecmditions 

(U.9)  5„(0)  =  . 

then 

(4ao) 

8„(x;e)  -  . 

are  solutions  of  t’roblen  B  with  f  *  P'  * 

nn 

To  invest l^te  the  solutions  of  (U*8)  we  employ  the 
analytleity  of  y,  s,  P  and  K  and  expand  them  in  the  format 
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(mn) 


(xje.S)  =  z(™>U;o)-!II  z'™’ (x;,^)  e''  , 

k=l  ^ 


(U.ll) 


T^  —  p  ^ j'  '  p ( V  —  e  \  f  ( nn )/,'\_k**l 

‘nin  ■  ^m  i ^  ^  ^  v  lr.)t 


""  ^  ^ 


where  we  have  used  the  last  of  (U.?)*  The  series  in  (^|•11) 
converge  in  soi-e  interval  lei  <  .  The  coefficients  in  the 

pji 

expansions  ai*e  deterriined  oy  substituting  (!^.ll)  into  (U.U) 
end  {L-,5J  •  Tt  linear  boundary  value  proLlens  for  y|^^^  and 

,  obtained  from  the  coefficients  of  e,  possess  solutions, 
if  and  only  if  the  inhonogenecus  terns  satisfy  the  appropriate 
orthogonality  condition.  This  condition  determines 
and  as 


(4*12) 


=  (-,2  .  2  A  .  -  A  ) 

1  mn  '  nn'  ^  inn  run  rvn^nm^  ’ 


p(nn),,,  ,  ('2,.2,-1(a_^,A„) 


whare. 


, - . .  3, c  c  ^ 

a  J^(f  )  '  nr'^mn  nn  nn  ^  nn  ^  ^nr.nn  * 


mn  2  I 


(U.X3) 


0 


2i 


Using  the  last  of  (i*.*!!),  the  bifurcation  equations 
reauc©  to 

(U.IU)  =  0  . 

If  CU.lii,)  have  continuous  solutions  ^  with 

^  IQ)  =  then/'®  Kuat  satisfy  the  bifurcation  conditions 
‘.'lan  inn  isn 

(Si.JS)  ='  ^  • 


Conversely,  if  satisfy  the  bifurcation  conditions  (U.l5) 
and  if 


(4*16) 


mn 


) 


d^' 


/  0 


ttien  ths  implicit  function  theorem  is  applicable  to  (^*11;) 

and  it  implies  that  the  bifurcation  equations  can  be  solved 

in  some  sufficiently  small  neighborhood  of  e  »  0. 

The  solutions  are  analytic  functi  is  for  suffi- 

ran 

ciently  small  c  and  stwisfy  6^„(0)  =  • 

ran  ran 

Inserting  i4»12)  and  (ii«13)  into  (1+.,15)  and  (4*18)  and 
defining  by 


(4.17) 


$ 


the  bifurcation  condi 'lions  reduce  to  the  cubles 
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(a.l8) 


^+E  ^  *Q  1  C+O  bO 

inn^  mn  an  3B3i  ^aa  ^  ' 


and  the  ’’solvability  conditions"  {i4.»l6}  reduce  to: 


(4.19) 


\  :< 


.  Aj 


^3 


1  _  \  3 


Here  we  have  used  the  notation. 


•;  ^4 


■■  ^%/V 


1 


nn 


2™|2-^-4a+  4)i  -3  -|-^  i  1  . 

'  nn  isffi  nn  y2/  :  ,,  ^nia[ 


Dsn  'ran' 


2(-.> 


(1^.^0) 


^am  ~ 


6 


(1+6^  ) 


flHil 


k 

¥ 


Thus  we  conclude  from  (4.18)  that  for  each  m  and  n  there  is 
at  least  one  real  5^  that  satisfies  the  bifurcation  condition 

4  ”  '^rm*  ^nin  roots  of  (4.l8). 

If  the  explicit  values  of  the  integrals  i  ,  defined 
in  (4*13)  ai*9  known  then  the  roots  of  (4.18)  cam  be  determined 
using  standard  formulae*  These  integrals  were  numerically 
determined  using  Simpson’s  Rule  with  500  mesh  points*  for  m 
and  n  in  the  range  1  ^  m,  n  <  25*  The  3essel  functions  were 
evaluated  using  the  procedure  of  [l9]  for  arguments  between 
zero  and  eight  and  a;;:propriate  asymptotic  formulae  were  used 

»  -  .  -  ,,,  __ 

Several  of  the  integrals  were  evaluated  using  1000  mesh 

points.  No  significant  difference  was  observed  with  the  re¬ 
sults  of  the  ^0  point  mesh.  Double  precision  arithmetic  was 
ea^Ioyed  and  all  calculations  were  preformed  on  the  IBM  7094 
computer  at  the  AEC  Computing  and  Applied  Mathematics  Center 
of  the  Courant  Institute  of  Mathematical  Sciences.  The  author 
is  Indebted  to  Dr,  P.  Bauer  for  conducting'  the  computation. 
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tor  larger  arguments.  V/e  found  that  ail  the  roots  of  (14., l8), 
for  m  and  n  in  the  ranged  1  <  m,  n  <  25»  ?»atisried  (It.  19) 
well  within  the  accuracy  of  the  computation.  In  fact,  the 
left  aide  of  (l4.«19)  was  usually  quite  large «  Hence  we  have 
shown  that  there  is  at  least  one  solution  of  Problem  B  for  P 
in  sorae  sufficiently  snail  neighborhood  of  each  double 
eigenvalue  with  m  and  n  in  the  abort  range,  Ve  conjecture 
that  there  are  solutions  near  each  double  eigenvalue  for  all 
K,n=l,2,,«,  •  We  also  find  that  for  the  following  indiciea 
(iul6)  h.is  three  real  roots: 


1) 

n^siafl. 

3<n<25»  n  odd 

2) 

n=w-3* 

lii<n<25»  n  odd 

3) 

ns=2»+l. 

2<n<2$, 

k) 

n»2B+3» 

65i<25, 

5) 

n=2HH-5, 

8<n<25* 

Thus  in  a  sufficiently  small  neighborliood  of  the  double 

eigenvalues  with  the  above  indicies  there  are  three  solutions 

of  Problem  B,  Hence,  as  ^ the  m-th  and  n-th  branches 

•  mn 

coalesce  and  when  P  »  p _  one  solution  nay  be  "destroyeo"  or 

‘  m 

a  third  solution  may  be  "created**  depending  on  the  values  of 
the  indiciea  m  and  n* 


T 


By  symetry  we  need  only  consider  n  >  m. 


2h 


5*  The  Xnterroediate  Buckling  Load 

It  is  convenient  to  refomulate  Probl‘“’Si  S  to  establish 
the  existence  of  the  interaediatw  buckling  load  cr 
equivalently  Xj^  =  ^?j.p  and  to  obtain  upper  and  lower  bounds 
on  Its  loagnitude*  F.quatiun  (2*lb)  is  integrated  and  the 
second  of  (2«3)  is  used  to  obtain 


X 

(5.1)  g*(x)  =  • 

'  0 

Then,  using  (2«Ii.b),  g(x)  is  given  as  a  functional  of  f(x)  by 


(5.2) 


- 


The  boundary  value  problem,  3,  is  reformulated  as  Problem  B* 
as  follows:  To  find  a  function  f(x)  which  possesses  a  contia* 
uous  second  derivative  and  satisfies  the  differential  equation 
(2.1a)  and  the  boundary  conditions 


(5.3)  f'(0)  =  f(l)  =  0  . 


The  function  g(x)  in  (2.1a)  is  defined  by  (5.1)  and  (5.2). 
These  equations  imply  that  g(x)  satisfies  (2.1b)  and 
g*(0)  =  g(l)  =  0.  This  statement  of  the  boundary  value 
problem  is  similar  to  the  one  introduced  in  [20l  in  a  study 
of  the  buckling  of  circular  plates. 

The  "energy"  functional  V  is  defined  as, 

1 

iS*h)  V[f(x)iX,f  ]  ay  [2(f *^(x)  -  Xf^(x))  ♦  g*^(x)]x^dx  , 

"0 
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whar^  g* (x)  ia  considered  as  a  functional  of  f{x)  defined  by 
(5*1) •  Thus  V  is  a  functicnai  of  f{x]  only  and  it  is  pro¬ 
portional  to  fcne  difference  between  the  potential  energies  of 
a  buckled  and  the  unbuckled  state. 

The  relationship  between  B*  and  V  is  obtained  by  first 
defining  a  class  of  functions.  f{x)  is  contained  in  A  or  is 
an  A-function  if  in  the  interval  0  ^  <  1  it  is  continuous 

and  satisfies  (5«3)»  f*(x)  is  an  function  and  all  integrals 
in  (5*1}«  (5*2)  and  (5*4/  exist*  Then  it  is  easy  to  show 
that  if  r{x)  ia  s  solution  of  B*  it  laakes  V  stationary  with 
respect  to  all  A-functions.  The  converse  of  this  result  can 
be  proved  using  the  methods  outlined  in  f20] •  The  converse 
states  that  if  f{x)  makes  V  stationary  (or  ndnimizes  V)  then 

/"(x)  is  continuous  and  f(x)  solves  B.  Here  f  (x)iA  is  said 

o 

to  ainisiize  V  for  a  fixed  X  and  p  if  V[f^(x)]  <  V[f  (x)]  for 
all  f(x) ^ A* 

The  result  conceding  the  existence  of  X^  is  now  stated 

as  the 


THEOREM.  If  for  every  finite  and  X  there  Is  an 
A-function  which  minimizes  V,  then  X^(r )  exists  and  it  is  in 

tha  intarral  )  <  X  ^  -P(-)  *  where  »(^ )  is  defined  in 

anj  i^  the  first  sero  of  Jj^(x)  *  0. 

The  theoreoi  is  a  direct  consequence  of  the  lennas  given 
below  and  the  following  inequality,  which  ia  easily  demon¬ 
strated  tieing  classical  methoda  [21] , 
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2 

(5.5)  /  frti)  -  >  0  ,  for  all  f(x)^  A  . 

0 

2 

Here  is  also  tb.e  lowest  eigenvalue  of  the  linearized 
buckling  problem  for  tbe  sy.netric  deformations  of  radially 
compressed  and  clamped  circular  plates: 

(5.6)  Gf{x)  >Xf(x)  =  0  ,  f*(C)  *  f(l)  s  0  . 


We  call  (5.6)  the  "equivalent"  circular  plate  prcblem.  The 
inequality  (5*5)  the  functional  V  given  in  (S.ii)  and  (5*1) 
immediately  yield 

LEI'iMA  1,  X  <  then  only  f  (x)  s  0  minimizes  V, 

The  form  of  the  functional  V  and  the  properties  of  the 
eigenvalues  (3.5a)  and  eigenvectors  (3*7)  of  the  linear  shell 
buckling  problem  yield 


LEIIMA  2.  If  X  =  ^(f)  then  there  exist  for  each  p  >  0 
a  function  f(x)e  A  such  that  v[f(x)]  <  0, 

The  proof  of  this  leraia  is  obtained  by  a  simple  calcula¬ 
tion  of  V’ff^(x)]  wnere  ^  eigenfunctions 

of  the  linear  theory  and  are  constants  restricted  to  lie 

n 

in  the  ranges  0  >  3  >  -2p  n  ■*  1,2,  •  The  constants 

n  n  n 

and  are  defined  by 


.2  ,  2 
— 


1  X 


J 

'0  [0 

I 


;  J  )x* dx^  x^^dx  , 
1  n 


1  2p- 


0  0 


jf(Ox’)x’dx' 
1  n 


X 

r 


j, (J  x’ )x’^dx’  x'^dx  >  0  . 
/  1  n 


27 


Lerma  1  iniplles  that  if  X  =  ^  then  g.l.b.  v[f{x)]  <  0.  The 
reiiflining  two  lecnaa  show  that  the  set  of  X’s  is  divided  into 
two  disjoint  sets:  those  \*s  for  which  the  mininiaQ  of  V  is 
zero  and  those  for  which  it  is  negative, 

L2JOiA  3-  IT  f(x)  =  0  nxini!n.ize3  V  for  X  =  X^,  then  onl? 
f (x)  5  0  minimizes  V  for  all  X  <  X^. 

Proof:  f  'ie  deduce  directly  from  the  form  of  7  that  if 
X  ^  X^,  V  >  0  for  all  A-functions.  Suppose  there  is  a 
X  »  X*  <  X^  and  an  A-functicn  f*(x)  4  0  which  minimizes  V, 
i«e«  V[f*(x);X*J  *  0.  Then  if  X  is  in  the  range  X'^  <  X  ^  X^, 
v[f^ixj  <  Yfr^iX*]  *  0.  This  is  in  contradiction  to  the  non- 
tkegatire  property  of  V  for  X  <  X^  and  the  proof  of  the  lemma 
la  eoMplete, 

LEMMA  1;.  If  the  minimum  of  7  is  negative  for  X  s  X^ 
than  it  is  negative_fQr  all  X  >  X^. 

Proof:  By  contradiction  using  Lemma  3» 

The  theorem  is  then  proved  if  we  define  X„(  ")  as  the 
l.u.b,  of  those  X  for  which  f(x)  -  C  mlnirizea  V.  Thus  for 

such  that  for  all  X  >  X^  there  arc 
buckled  states  with  less  energy  than  the  unbuckled  state  and 
for  X  <  X^  all  buckled  states  have  greater  energy  than  the 
tmbuckled  state. 

TtM  intemediate  load  can  also  be  characterized  by  the 
^  This  form  of  the  proof  was  suggested  by  Dr.  M.  Mevnai^ 


each  there  is  e  X  «  X^ 


I 
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miniaiuiR  property, 


(5.7) 

^  f(x)  ^ 


(2/*^+  g‘^}x^dx 


0  t 


0 


2f^x^dLX 


where  g*(x)  is  d-fiaed  ia  (5.1}.  'Joper  bounds  for  are 
obtained  by  selecting  trial  A-f’-mct ions  to  r.ake  tt*e  quotient 
in  {5.7)  as  snail  as  possible.  For  example,  functions  of  tine 
form  f(x)  =  3F(x}  are  considered,  where  P(x)  -  A  is  a  specified 
function  and  3  is  an  arbitrary  constant.  The  quotient  is 
then  a  function  of  5  only  and  -3  is  determined  to  minimize  it. 
Some  of  the  trial  functions  that  were  used  and  tne  resulting 
uDpcr  bounds,  X,,  A  )  =  a,  are  shown*  in  Table  I. 

^  s w X 


;  i 

t  ; 

■  i 

i  ' 

1 

:  -  V- 

.  f  * 

1 

i  i  ^ 

] 

?{x) 

*1 

h 

r~ 

:  I 

l-x  1 

15-C 

b.5io6 

2 

f  « 

i'-X 

Ic  .C 

3.7037 

3 

3 

I*x-^ 

1  *’■  Cp 

2  - 

k 

1-x^ 

i  .  20 

1.7166 

5 

2 

1+-X-2X 

16,7576 

’  2.7636 

6 

2(2-x-x^) 

1  u.cOSS 

The  calculations  were  perfornect  by  Szeto. 
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The  results  are  suBzsarized  in  Pig.  3  where  we  also  show 
accurate  approximations  of  obtained  frcan  a  numerical 
solution  of  Problem  B.  Details  of  these  and  other  calcula* 
tions  will  be  reported  elsewhere  [22] . 

6«  Other  Bifurcation  Problems 

Analogous  results  can  be  obtained  for  other  bifurcaticn 
buckling  problems.  For  exaiaple,  for  the  boundary  valve 
problem  consisting  of  (2«1)9  (2.3)  (2.5)  which  we  call 

Problem  3^>  we  can  show  that 

where  X  is  now  the  rninimum  eigenvalite  of  the  linearized  shell 
buckling  theory  using  the  boundary  conditions  (2.5)  in  place 
of  (2.1i.).  Upper  bounds  for  ^(p)  ai*e  obtained  from  a  formula 
similar  to  (5*7) .  These  results  are  in  Pig.  3  'rlth 

the  predictions  of  X^i  obtained  fron  a  numerical  solution 
of  the  boundary  \^lue  problem  [22]. 
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Figiir©  1 

Plgui*©  Z 

Figure  3 


Captlona  for  Figure 3 

Sketch  of  conjecti^red  load  defonoation  curve  for 
bifurcation  buckling.  Here  D  is  a  representative 
def 'vrffiatlcn,  e.g.  tne  normal  displacesKsut  of  Ltae 
cap^s  center* 

Shell  geometry. 

Intemediate  buckling  loads  for  Problems  B  and  B*» 
Tha  upper  bound  curves  are  essentially  the  envelopes 
of  the  curves  given  in  Tables  I  and  II, 
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